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Abstract: We report on a measurement of the neutron deteetion effieieney in Nal crystals in the 
Crystal Ball detector obtained from a study of single TtP photoproduction on deuterium using the 
tagged photon beam at the Mainz Microtron. The results were obtained up to a neutron energy of 
400 MeV. They are compared to previous measurements made more than 15 years ago at the pion 
beam at the BNL AGS. 
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1. Introduction 

The A2 Collaboration at MAMI is engaged in a program ||I]] to study yn-interactions for which 
neutron detection is of crucial importance. The A2 experiment consists of several detector systems 
described below. Two of them, the Crystal Ball (CB) and TAPS, have a good neutron detection ef¬ 
ficiency. The Nal CB detector was built originally at SLAC in the mid 1970s for use with colliding 
e+e^ beams at SPEAR |^. Later the CB was involved in experiments at DORIS and Brookhaven 
National Laboratory before arriving at MAMI where it demonstrates an excellent performance for 
photon and charged particle detection. Neutron detection with this detector has been studied less. 
The CB neutron detection efficiency was measured previously in a 1997-1998 run at BNL using 
the reaction n^p ^ K^n. There are very few references in the literature to neutron interactions in 
Nal. A short summary of this topic can be found in [Q]. Previous measurements are not directly ap¬ 
plicable to the current status of the CB because of possible aging effects, differences in thresholds 
and in methods of analyzing the data. In this paper, we describe a new study of the CB neutron 
detection efficiency based on measuremenfs of single 71° photoproduction on a liquid deuterium 
target using the tagged photon beam at MAMI. 
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Figure 1. Layout of the A2 detector systems. 


2. Experimental setup 

The measurements were performed at the tagged photon facility of the Mainz Microtron acceler¬ 
ator (MAMI) [Q]. An electron beam was used to produce bremsstrahlung photons, which were 
tagged with the upgraded Glasgow magnetic spectrometer ||^. The target was a Kapton cylinder 
of 4 cm diameter and 10 cm length filled with liquid deuterium. Photons, charged pions, and re¬ 
coil nucleons, produced in the target, were detected with an almost 4?! electromagnetic calorimeter 
schematically shown in figure |I|. If combined fhe CB and fhe forward TAPS defector |^, alfhough 
fhe laffer was nof used in fhe currenf analysis. 

The CB defecfor is a sphere consisfing of 672 opfically isolated Nal(Tl) crysfals, shaped as 
fruncafed friangular pyramids, which poinf foward fhe cenfer of fhe sphere. The crysfals are ar¬ 
ranged in fwo hemispheres fhaf cover 93% of An sr or polar angles of 20-160°, sifting oufside a 
cenfral spherical cavify wifh a radius of 25 cm. This cavify holds fhe fargef and inner defectors. 
The Nal crysfals have a lengfh of 40.7 cm, which is equal fo 15.7 radiation lengfhs or ~1 hadron 
inferacfion lengfh. The fargef was locafed in fhe cenfer of fhe CB and was surrounded by a Par- 
ficle Idenfificafion Defecfor (PID) used to disfinguish befween neufral and differenl charged 
particles based on dE/dx measuremenfs. If was made of 24 scintillator bars (50 cm long, 4 mm 
fhick) arranged as a horizonfal cylinder wifh a radius of 6 cm. The PID was surrounded by fwo 
Cylindrical Multiwire Proporfional Chambers (MWPC). The inner (outer) chamber has a radius of 
7.4 (9.45) cm and a lengfh of 57 cm. Each chamber measured fhe fhree-dimensional coordinates of 
a charged-parficle frack as a resulf of a readouf of fhree signals, one from a horizonfal anode sense 
wire and fhe ofher fwo from spiral sfrip cafhodes 0. 
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3. Overview of the method 


For the neutron detection efficiency measurement, we chose the 71** photoproduction on deuterium, 

y+d ^ Tip+ p + n, (3.1) 

where the tP decays to two photons. The reaction kinematics is completely determined if the 
pion and proton are detected and the beam energy is known. The tP momentum is reconstructed 
using the energies and directions of the decay photons. The proton momentum is determined by 
its deposited energy in the CB and track direction in the MWPC and the photon beam energy is 
given by the tagging spectrometer. The reaction vertex was determined from the intersection of the 
proton track with the photon beam axis and allowed rejection of events from target walls. Then 
reaction ( |3.1D can be identified by selecfing fhe neufron by missing mass. This procedure rejecfs 
evenfs with more than one final-state pion and provides a reconstruction of the neutron momentum 
vector. If neutron momentum points to the CB (polar angles 30 - 150°), then a neutral hit in the 
CB in the same direction is searched for. Then the neutron detection efficiency is calculafed as fhe 
rafio of regisfered neufral hifs fo all neufrons. The indicafed angular range is chosen fo be smaller 
fhan the full acceptance of the CB to eliminate edge regions that inevitably have efficiency losses 
due fo uncerfainfies in predicfed and real hif positions. 

Aparf from being kinematically complefely defermined, reaction ( ^T| ) has other advantages. It 
has a large cross section and a broad energy and angle spectrum of neutrons. The energy spectrum 
extends from nearly zero, when the pion is produced on a proton with a spectator neutron, to the 
maximum energy from pion production on the neutron. The latter was limited by the beam energy 
and the available statistics. Reaction ( ^d] ) was identified from four clusfers in fhe CB and one hif 
in fhe PID. Two clusfers arise from tP decay phofons, fhe fhird one from proton energy loss, and 
fhe fourlh from fhe neufron inferacfion wifh Nal. A clusfer is defined as a group of nof more fhan 
13 adjacenf Nal crysfals each wifh deposifed energy larger fhan 2 MeV centered around fhe crysfal 
wifh maximum energy deposifion. The minimum cluster summed energy was sef fo 15 MeV in fhe 
analysis. In fhis measuremenf, we selected evenfs wifh fhree or four clusfers in the CB and only 
one hit in PID to reduce background from multiple pion production. In the following sections, we 
describe all above mentioned steps in more detail. 

4. Photon Tagger 

For the neutron detection efficiency measuremenf, we used dafa from fhe deuferium run faken af 
MAMI [Q] in March 2013. Bremssfrahlung phofons, produced by fhe 883 MeV elecfrons in a 10 
pm Cu foil and collimated by a Pb collimator 3 mm in diamefer, were incidenf on a 10 cm long 
liquid deuferium fargef locafed af the center of the CB. The beam was 1 cm in diameter at the target. 
The energies of the incident photons were measured by detecting the post-bremsstrahlung electrons 
in the 352 channel focal-plane detector of the Glasgow-Mainz tagger [^. The energy resolution of 
the tagged photons is determined mostly by the width of the tagger focal-plane detectors and by 
the electron-beam energy. For the present beam energies, the typical width of a tagger channel was 
about 2 MeV. The data were taken with a trigger that required the total energy deposit in the CB 
to exceed 40 MeV. Figure ^ shows the time distribution of tagged photons relative to the trigger. 
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Figure 2. Time distribution of tagged photons. Figure 3. Invariant 77-mass. A Gaussian fit to the 

peak of the distribution is shown in red. 

The peak near 0 ns is due to the prompt photons that produced the triggers. The entries outside 
and under this peak are due to random photons. The r.m.s. width of this peak is equal to 2.7 ns 
and is determined mostly by the intrinsic time resolution of the CB and time alignment of different 
channels. More accurate adjustment of the channel time offsets can improve the time resolution 
by a factor of two but it was not important for this analysis. We selected the prompt photons in 
a time window (-10 ns, -rlO ns) and subtracted background produced by random photons for all 
distributions where the photon energy was used. The wide time range accepted by the tagger made 
it possible to use wide time windows for random photons (-230 ns to -30 ns and 30 ns to 230 ns) 
that resulted in a negligible subtraction error. 

5. Selection of the single neutral pion 

The first step for reconstruction of reaction ( pd] ) was neutral pion selection. This was done by 
choosing the photon pair, measured in the CB, with an invariant mass niyy closest to the 7i^ mass 
nijfi. An additional cut on the sum of energies of the two photons Eyy > rejects low-energy 
photons. The tuyy distribution obtained with this cut is shown in figure ^ wifh a gaussian fif fo fhe 
peak. The disfribufion demonsfrafes fhe excellenf resolufion of fhe Nal CB defector in m^o-mass 
of 9 MeV (r.m.s.). For tt** selection, we used a 120 < niyy < 150 MeV cuf. Identification of fhe 
proton and measuremenf of ifs energy were done in a few sfeps. The prof on frack was required fo 
have hifs in fhe PID and fhe CB and space poinfs in fwo MWPC compatible wifh a sfraighf line 
fhaf has a poinf of closesl approach less fhan 10 mm from fhe beam axis wifhin fhe volume of fhe 
fargef. The separation of protons from pions and elecfrons/posifrons is displayed in figure ^ where 
fhe CB clusfer energy is ploffed along fhe horizonfal axis and dE/dx is ploffed along fhe verfical 
axis. dE/dx = FpiosinO, where FpiD is fhe energy deposifed in fhe PID and Q is fhe frack angle 
wifh respecf fo fhe beam axis. The facfor sin Q correcfs fhe deposifed energy for frack lengfh in 
fhe PID scinfillafor segmenfs which are parallel fo fhe beam axis. lonizafion losses for profons 
depend on fheir energies. In figure profons lie above fhe empirically drawn line while pions and 
elecfrons/posifrons lie below. 
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Figure 4. Pion/proton separation on the two dimen¬ 
sional plot dE/dx as measured by PID vs. CB to¬ 
tal energy. Protons lie above the empirically drawn 
line, pions - below. 



CB Energy, MeV 

Figure 5. Proton energy correction for ionization 
losses and Nal light yield. 7p - proton kinetic energy 
at interaction point, E - energy measured in the CB 
detector, 0 - polar angle of proton track. Line - htted 
correction curve. 


Two effects must be taken into account to determine the proton energy at the interaction point. 
First, due to higher dE/dx for protons that are stopped in the Nal crystals of the CB detector, the 
energy calibration for protons is different from that of minimum-ionizing electromagnetic showers. 
Second, the measured proton energy in the CB has to be corrected for ionization losses in detector 
materials between the interaction point and the CB crystals. These materials are 2 cm of liquid 
deuterium, 1 mm CH of the target walls, 4 mm PID scintillator, and 1.5 mm stainless steel of 
the CB spherical inner wall. The proton energy at the interaction point was found by using the 
correction curve shown in figure ^ where (Tp — £')sin0 is plotted vs. the energy {E) measured 
in the CB. Here Tp is the proton kinetic energy at the interaction point calculated from energy 
and momentum conservation for reaction ( p.l| ) using photon and pion momentum vectors and the 
direction of the proton track. At each bin of CB energy, (7p — £')sin0 was fitted with a gaussian 
and the maximum was found. Then these maxima at different energies were fitted with the function 
(exp(ci — C 2 E) — C 3 ) with three free parameters c,, / = 1 — 3. The efficient performance of this 
procedure is demonstrated in figure ^ where distributions of the events over missing mass 
in reaction ( ^ ) are given for a few intervals of the beam energy. Missing mass is defined as 
^miss ~ \l (-^r + -^d “ ~ -^p)^’ where Pi (/ = 7 ,d, K^,p) are the four-momentum vectors for the 

incident photon, target deuteron, neutral pion, and proton, respectively. Good agreement between 
the peak positions and the neutron mass is clearly seen in all four photon energy ranges. 

6. CB neutron detection efficiency 

Neutrons were selected by a cut on the missing mass M^iss centered at the mass of a neutron 
890 < Mmiss < 990 MeV. The momentum vector for a selected neutron was calculated using energy 
and momentum conservation. This vector must point to a CB region slightly smaller than instru¬ 
mented with Nal crystals to eliminate edge effects; that is, the angle Q of this vector with respect 
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Figure 6 . Missing mass distributions for the reaction d ( 7 , n^p)X for different incident photon energies. The 
peak position demonstrates good agreement with the neutron mass at all photon energies. 


to the beam direction must be 30° < 6 < 150°. The angles A between this vector and the direction 
vectors for neutral hits in the CB were determined. Hits from protons and photons from 71 ^ decay 
were excluded from this procedure. The distribution over cos A is given in figure ^ It has a large 
peak at A = 0 signifying the real neutron hit in the Nal crystals and a small background of random 
hits. We chose cos A = 0.85 as the boundary between the region of real neutron hits in the CB and 
background. The small uncertainty connected with this cut will be discussed later. 

The energy distribution of neutrons from reaction ( |3.lD identified by missing mass and incidenf 
on fhe CB is shown by fhe solid line in figure p. Those neutrons that gave hits in the CB with cos A 
< 0.85 have the energy distribution shown by the dashed line. The first distribution is peaked at low 
energy, which demonstrates the large contribution of spectator neutrons from pion photoproduction 
on the proton. The second distribution goes to zero at low energies because of the low detection 
efficiency for neufrons. The ratio of fhese fwo disfributions gives fhe CB defecfion efficiency for 
neufrons shown in figure ^ and lisfed in Table |T] for 15 and 20 MeV clusfer energy fhresholds. The 
defection efficiency is averaged over ±10 MeV relative fo fhe lisfed energy values. The sfafisfical 
and fhe fwo main sysfemafic uncerfainfies for 15 MeV fhreshold are also given. The sfafistical 
uncerfainfies vary from less fhan 0.1% af fhe lowesf energy fo 3% af 400 MeV. The increase resulfs 
from fhe decrease in cross section for reacfion (^) as proton and neufron energy increases. The 
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Figure 7. Cosine of the angle between the 
predicted neutron momentum vector and the 
direction vector for third neutral hit in CB. 



Figure 8. Kinetic energy for neutrons identi¬ 
fied by missing mass - solid histogram. When 
a neutron hit is registered the dashed his¬ 
togram is obtained. 



Figure 9. CB detection efficiency for neutrons as a function of neutron kinetic energy. Downward (upward) 
triangles - present measurements for a 15(20) MeV cluster energy threshold; circles - previous measurement 


at BNL in 1997-98 for a 20 MeV threshold [|]. 


systematic uncertainty I was estimated by changing the cut on missing mass from (890 - 990) MeV 
to (890 - 1010) MeV. The wider window gives slightly lower detection efficiency due to a small 
contamination from two-pion production. The smaller value of the lower boundary of the cut has 
negligible effect on the neutron detection efficiency. The sysfemafic uncerfainfy II was esfimafed 
by changing the cut on cos A from the standard one cos A > 0.85 to cos A > 0.65. Changing this 
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Table 1. Neutron detection efficiency as a function of neutron kinetic energy. Column 2 gives the efficiency 
and statistical uncertainty for a 20 MeV threshold, column 3 gives the efficiency and statistical uncertainty 
for a 15 MeV threshold, columns 4 and 5 give systematic uncertainties I and II, and column 6 gives the total 
systematic uncertainty. See text for details. 


Energy, 

MeV 

Efficiency, % 

(20 MeV threshold) 

Efficiency, % 

(15 MeV threshold) 

Sys. 

I, % 

Sys. 

n, % 

Sys. 
total, % 

20 

0.69±0.01 

1.39±0.01 

0.05 

0.02 

0.05 

40 

6 .6±0.1 

9.7±0.1 

0.1 

0.1 

0.2 

60 

14.4±0.1 

18.8±0.1 

0.1 

0.2 

0.2 

80 

20 .6±0.2 

25.3±0.2 

0.1 

0.3 

0.3 

100 

25.7±0.2 

30.3±0.2 

0.1 

0.3 

0.3 

120 

29.0±0.3 

33.1±0.3 

0.2 

0.3 

0.3 

140 

31.7±0.3 

35.2±0.3 

0.4 

0.4 

0.5 

160 

33.5±0.4 

36.6±0.4 

0.4 

0.4 

0.6 

180 

35.2±0.5 

38.1±0.5 

0.5 

0.4 

0.7 

200 

36.1±0.5 

38.5±0.5 

0.5 

0.5 

0.7 

220 

36.6±0.6 

38.8±0.6 

0.7 

0.4 

0.8 

240 

36.8±0.7 

38.5±0.7 

0.6 

0.4 

0.7 

260 

38.9±0.9 

40.8±0.9 

0.4 

0.4 

0.6 

280 

37.8±1.0 

39.5±1.0 

0.7 

0.7 

1.0 

300 

38.1±1.1 

40.1±1.1 

0.8 

0.5 

1.0 

320 

39.1±1.4 

41.0±1.2 

0.6 

0.7 

1.0 

340 

36.8±1.7 

38.6±1.7 

0.8 

0.6 

1.0 

360 

41.1±2.1 

42.0±2.1 

0.5 

0.4 

0.6 

380 

38.8±2.6 

40.3±2.6 

0.7 

0.9 

1.2 

400 

39.8±3.4 

40.6±3.4 

0.4 

0.3 

0.7 


cut resulted in a higher detection efficiency as expected from the larger acceptance for neutron 
hits in the CB. Other cuts were also tested: tt** - mass, window for prompt photons, position of 
interaction point in the liquid deuterium target, exclusion of dead channels (1% of all channels). 
Each of these cuts had an order-of-magnitude smaller effect on neutron detection efficiency than 
I and II, mentioned above. The total systematic uncertainty was obtained as the quadrature sum 
of systematic uncertainties I and II and it is given in the last column of Table [I[ For the most 
part, its value is close to the statistical uncertainty. Table [I] also demonstrates the sensitivity of 
neutron detection efficiency to the cluster energy threshold. An increase of the threshold from 15 
to 20 MeV results in a relative decrease of neutron detection efficiency by 30% at 40 MeV neutron 
energy and by less than 5% for the energies higher than 300 MeV. In figure the CB neutron 
detection efficiency measured more than 15 years ago [Q] when the detector was in operation at 
BNL is also shown. The shape of its energy dependence for a 20 MeV threshold is in reasonable 
agreement with the present measurement for a 15 MeV threshold. Absolute values of CB neutron 
detection efficiency for the present measurement with a 20 MeV threshold are 2-4% smaller than 
measured at BNL. 











Figure 10. Energy deposited in the CB for three in¬ 
cident neutron energies (marked by arrows): 350 ± 
30 MeV - solid black line, 200 ± 20 MeV - dashed 
red line, 50 ± 10 MeV - dash-dot blue line. The 
histograms were multiplied by the given factors to 
compensate for the strong energy dependence of the 
neutron flux. 


Figure 11. Energy weighted “radius" of the clus¬ 
ter for proton (black points), neutron (red upward 
triangles) and photon (blue downward triangles) in¬ 
duced showers. The peak at zero is due to clusters 
that consist of only one crystal. 


7. Neutron interaction characteristics 


Neutrons interacting in the CB, which is not designed to be a totally absorbing hadronic calorimeter, 
generally deposit only part of their kinetic energy. This fact is illustrated by figure |I^ where 
distributions over cluster energies recorded in the CB for hits identified as neufrons are given for 
fhree selecfed infervals of incidenf neufron energy: 50 ± 10 MeV, 200 ± 20 MeV, and 350 ± 
30 MeV. For each incidenf energy, fhe response in fhe CB is disfribufed from fhreshold fo a maximal 
energy fhaf slighfly exceeds fhe incidenf energy. This effecf can be affribufed fo energy resolufion 
and fo some difference in calibration befween hadron-induced and phofon-induced showers. The 
laffer was used fo defermine fhe neufron response. Figure E shows fhaf fhe energy in a clusfer 
induced by a neufron has limifed ufilify for an esfimafion of fhe neufron energy; however, fhe shape 
of a clusfer can be more informafive. In figure normalized disfribufions over clusfer “radius" 
are given. The clusfer “radius" is an energy-weighfed sum of disfances befween fhe crysfal wifh 
maximum energy and each remaining crysfal in a clusfer. The clusfer “radius" is equal fo zero for 
a clusfer composed of only one crysfal. For a clusfer induced by a neufron fhe probabilify fo have 
zero “radius" is nearly 30%. If is much higher fhan for a phofon clusfer, for which fhe probabilify 
is only 4%. This difference can be used fo discriminate befween neufron- and phofon-induced 
showers. Profon-induced showers have an even larger probabilify for zero “radius" fhan neufron- 
induced showers. This information can also be used fo confrol a sample of profons. 


8. Discussion and conclusion 

We have measured fhe defecfion efficiency in fhe Nal CB defector for neufrons in fhe kinetic en- 
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ergy range 20 - 400 MeV using 71° photoproduction on deuterium. A previous measurement using 
the same detector was performed at BNL ||3|] in 1997-1998 with the 71^p n^n reaction. The 
present measurements agree closely with the previous measurements for the shape of the energy 
dependence of neutron detection efficiency but we find efficiencies fhaf are 2-4% smaller. This 
effecf cannof be unambiguously explained as a resulf of possible degradation of fhe lighf yield of 
Nal crysfals, which was pari of fhe mofivalion for Ihis sludy. Many hardware and soflware changes 
have occurred belween Ihese Iwo measuremenls. Differenl energy calibration and cluslering algo- 
rilhms can also resulf in a slighfly differenl neufron defecfion efficiency. The presenf measuremenf 
gives fhe up-lo-dale sfalus of neufron defection wilh fhe CB delecfor and can be used in various 
physics analyses fhaf involve neulrons and for improving fhe description of neufron defecfion in 
Monle Carlo simulafions. 
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